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ABSTRACT 



We investigate the EUV and X-ray flare rate distribution in radiated energy of 
the late-type active star AD Leo. Occurrence rates of solar flares have previously 
been found to be distributed in energy according to a power law, dN/dE oc E~ a , 
with a power-law index a in the range 1.5—2.6. If a > 2, then an extrapolation of 
the flare distribution to low flare energies may be sufficient to heat the complete 
observable X-ray/EUV corona. 
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We have obtained long observations of AD Leo with the EUVE and Bep- 
poSAX satellites. Numerous flares have been detected, ranging over almost two 
orders of magnitude in their radiated energy. We compare the observed light 
curves with light curves synthesized from model flares that are distributed in 
energy according to a power law with selectable index a. Two methods are ap- 
plied, the first comparing flux distributions of the binned data, and the second 
using the distributions of photon arrival time differences in the unbinned data 
(for EUVE). Subsets of the light curves are tested individually, and the quiescent 
flux has optionally been treated as a superposition of flares from the same flare 
distribution. We find acceptable a values between 2.0—2.5 for the EUVE DS 
and the BeppoSAX LECS data. Some variation is found depending on whether 
or not a strong and long-lasting flare occurring in the EUVE data is included. 
The BeppoSAX MECS data indicate a somewhat shallower energy distribution 
(smaller a) than the simultaneously observed LECS data, which is attributed to 
the harder range of sensitivity of the MECS detector and the increasing peak 
temperatures of flares with increasing total (radiative) energy. The results sug- 
gest that flares can play an important role in the energy release of this active 
corona. We discuss caveats related to time variability, total energy, and multiple 
power-law distributions. Studying the limiting case of a corona that is entirely 
heated by a population of flares, we derive an expression for the time-averaged 
coronal differential emission measure distribution (DEM) that can be used as a 
diagnostic for the flare energy distribution. The shape of the analytical DEM 
agrees with previously published DEMs from observations of active stars. 

Subject headings: Stars: activity — stars: coronae — stars: flare — stars: individual 
(AD Leo) — stars: late-type 

1. Introduction 

The physics of coronal heating remains one of the most fundamental problems in stellar 
(and solar) astrophysics. The subject has been reviewed extensively from the point of view of 
theoretical concepts (Ionson 1985; Narain & Ulmschneider 1990; Zirker 1993), observational 
solar physics (e.g., Benz 1994), and stellar physics (e.g., Haisch & Schmitt 1996), where the 
cited work stands exemplary for a large body of literature available. It is somewhat surprising 
that the nature of the "coronal heating mechanism(s)" still eludes agreement given high- 
resolution imaging of solar coronal structures or large statistical samples of stellar coronal 
X-ray observations. For example, there is no unequivocal agreement on whether all, or any, 



- 3- 



of the X-ray coronal energy detected from certain classes of stars is magnetic in origin. 

Coronal heating is of particular interest to stellar astrophysics since it relates directly to 
our understanding of coronal structure and dynamics, information that is usually obtained 
by means of indirect modeling. Apart from heating models involving acoustic heating, for 
example on F-type stars (see Mullan & Cheng 1994; although the resulting X-ray flux would 
be much smaller than observed - see St§pieh & Ulmschneider 1989), the currently advocated 
mechanisms are of two types: (i) Steady heating mechanisms, e.g., by steady electric current 
dissipation or MHD waves, and (ii) heating by explosive energy release, e.g., coronal flares. 
The latter are attractive heating agents since flares do heat plasma efficiently, although only 
episodically since the radiative and conductive losses rapidly cool plasma to pre-flare levels, 
typically within minutes to hours. 

The flare heating hypothesis has gained momentum in particular from solar, but also 
from stellar observations during recent years. If the quasi-steady ( "quiescent" ) coronal emis- 
sion is to be explained by flare contributions, flares must act as stochastic heating agents. 
Parker (1988) proposed that shuffling of magnetic field footpoints in the photosphere by the 
convective motions leads to tangled magnetic field lines in the corona and thus to current 
sheets. With increasing winding of magnetic fields, the necessary energy may be transported 
into the coronal magnetic field where it is released by sudden relaxation involving reconnec- 
tion. Parker estimates that energy dissipation occurs in packets involving 10 24 — 10 25 ergs 
( "nanoflares" ) . The flare-heating hypothesis resolves to the basic question of whether or 
not the statistical ensemble of flares (in time and energy) suffices to heat the apparently 
nonflaring coronae. 



2. Statistical Flare Observations: A Brief Overview 

Solar observations have provided evidence for the presence of numerous small-scale flare 
events occurring in the solar corona at any time (Lin et al. 1984; Porter, Fontenla, & Simnett 
1995; Gary, Hartl, & Shimizu 1997; Krucker & Benz 1998). Their distribution in energy was 
found to be a power law (Datlowe, Elcan, & Hudson 1974; Lin et al. 1984) of the form 

where dN is the number of flares (per unit time) with a total energy (thermal or radiated) 
in the interval [E,E + alE\. The power-law index a is crucial: If a > 2, then the energy 
integration (for a given time interval), 
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(assuming E min <C E max and a > 2 for the last approximation) diverges for E min — > 0, i.e., 
by extrapolating the power law to sufficiently small flare energies (microflares with ~ 10 27 — 
10 30 erg, nanoflares with ~ 10 24 — 10 27 erg), any energy release power can be attained (Hudson 
1991). This is not the case for a < 2. Evidently, then, one needs to measure the energy 
distribution of a statistically relevant number of flares. Solar studies have repeatedly resulted 
in a values on the order of 1.6-1.8 for ordinary solar flares (Crosby, Aschwanden, & Dennis 
1993), insufficient for a coronal microflare- heating model. On the other hand, such values 
have proven interesting for statistical modeling since they are also found from avalanche 
models (Lu & Hamilton 1991). However, the implicit assumption is that the same power 
law continues to very low energies with unchanged a. If a steepens for low energies, then 
the microflare-heating concept may well be in order (Hudson 1991), and such a steepening is 
supported by simulations of avalanche models for small events (Vlahos et al. 1995). Detailed 
studies of small events in the soft X-ray (SXR) and extreme ultraviolet (EUV) ranges with 
Yohkoh, SoHO, and TRACE still provide a somewhat ambiguous picture (Shimizu 1995; 
Krucker & Benz 1998; Aschwanden et al. 2000; Parnell & Jupp 2000), although recent 
statistical investigations suggest a = 2.0 — 2.6 for small flares in the quiet solar corona 
(Krucker & Benz 1998; Parnell & Jupp 2000). An extrapolation to unseen microevents is 
always required to explain all of the coronal energy release, in some cases by several orders 
of magnitude. Aschwanden et al. (2000) argued that the extrapolation is limited if small 
flares follow some scaling laws established at higher energies, i.e., their temperatures drop to 
sub-coronal values, their densities become unreasonably small, and their surface area filling 
factor approaches unity possibly above the low-energy limit required. 

There are two critical questions to ask here: (i) Is the power-law index of the flare 
energy distribution different for different flare energy ranges? Even if a < 2 for ordinary 
flares, a steepening at lower energies to a > 2 can have important consequences for the 
coronal heating energy budget, and the recent findings for microflares seem to support this 
hypothesis, (ii) Does the derived power-law index a depend on the spectral range in which 
the observations were made? Incidentally, the recent microflare observations that result in 
a > 2 were mostly obtained in the EUV range where the dominant emission comes from hot 
plasma. It is by no means evident that the hard X-ray production is proportional to the 
thermal flare energy (Feldman, Doschek, & Klimchuk 1997), and it is clear that the energy 
released in the hard X-ray (HXR) range (that has often been used for flare statistics) is a 
very small fraction of the total flare energy. Heating and particle acceleration systematics 
in the flare energy release may affect both problem areas: The flare temperature and the 
production of non-thermal particles systematically depend on the total flare energy (e.g., 
Feldman, Laming, & Doschek 1995; Porter et al. 1995; Aschwanden 1999; Krucker & Benz 
2000). We will address both questions in the course of the present study. A summary of a 
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few relevant observations of solar flare energy distributions is given in Table 1. A further 
important problem visible in this table relates to peak flux measurements: These are reliable 
indicators for the total flare energy only if the flare duration does not depend on the total 
flare energy (Crosby et al. 1993; Porter et al. 1995). 

There has been some interesting progress also on the stellar side. Most of the evidence for 
an important role of flares relates to magnetically active stars with coronal energy loss rates 
orders of magnitude larger than the Sun's. Doyle & Butler (1985) and Skumanich (1985) 
noted that the quiescent coronal X-ray luminosity of active stars, Lx, is correlated with their 
time-averaged U-band flare luminosity, which may indicate that the statistical ensemble of 
flares is responsible for the energy deposition in the corona. Analogously, Audard et al. 
(2000) find that the X-ray flare frequency (above a given energy limit) also correlates with 
L x . Wood et al. (1996) interpret broadened ultraviolet-line wings as being due to a large 
number of transient, explosive events similar to those previously identified in the solar corona 
(Dere, Bartoe, & Brueckner 1989). Observations of active stars often show slow variability 
within a factor of 2 over both short and long terms, apart from obvious flaring (Ambruster, 
Sciortino, & Golub 1987; Collura, Pasquini, & Schmitt 1988; Giidel et al. 1996; Kashyap & 
Drake 1999; Sciortino et al. 1999). Statistically, the hotter X-ray emitting plasma component 
is more variable than the cooler one, indicating that flares may be involved (Giampapa et al. 
1996). Further, small flare events with energies of the order of 10 27 — 10 28 ergs have become 
observable with the Hubble Space Telescope in cool M dwarfs, and some of their energy 
distributions suggest a > 2 (Robinson et al. 1995, 1999, 2001) although we note that these 
transition region events may not be related to the physics of coronal heating. 

Studies of statistical flare energy distributions in stellar coronae have been rare, due to 
the paucity of relevant data sets (Table 1). Collura et al. (1988) and Pallavicini, Tagliaferri, 
& Stella (1990) derived a = 1.5 — 1.7 for a sample of M-dwarf flare observations with 
EXOSAT, and Osten & Brown (1999) report a = 1.6 for a set of flares on RS CVn binary 
systems observed with EUVE. All observations obtained flare energies simply by integrating 
count rates during flares selected from the light curves. Also, these investigations collected 
flare information from numerous different stars, regardless of distance (introducing selection 
bias for small flares), luminosity (and hence overall magnetic activity level), binarity, etc. 
To avoid possible selection bias, Audard, Giidel, & Guinan (1999) and Audard et al. (2000) 
applied a flare search algorithm to EUVE light curves of individual active main-sequence 
stars, taking into account flare superpositions and various binning to recognize weak flares. 
Although the statistical confidence intervals are broad, the results indicate a predominance 
of relatively steep power laws including a > 2, suggesting that the extension of the flare 
distribution to lower energies may contribute considerably to the detected energy loss. 
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Table 1. PREVIOUS MEASUREMENTS OF FLARE ENERGY DISTRIBUTIONS 

(SELECTION) 



Measurement 


Photon energy range 


Flare energy 


a 


Reference 


O (jflQj I m 










Thermal energy 


EUV Fc IX, X, XII 


10 25 - 3 X 10 26 erg 


2.3-2.6 


Krucker & Benz 1998 


Thermal energy 


EUV Fc IX, X, XII 


3 X 10 23 - 10 26 erg 


2.0-2.6 


Parnell & Jupp 2000 


Thermal energy 


EUV Fe IX, X, XII 


10 24 - 2 X 10 26 erg 


1.79 ±0.08 


Aschwanden et al. 2000 


Thermal and radiated energy 


SXR 


lO 2 ? - 10 29 erg 


1.5-1.6 


Shimizu 1995 


Peak flux 


SXR 


10 25 - 10 29 erg 


1.4-1.8 


Shimizu 1995 


Peak flux 


1-6 kcV 


normal flares 


1.75 


Drake 1971 


Radiated energy 


1-6 keV 


normal flares 


1.44 


Drake 1971 


Peak flux 


1.5-12 keV 


GOES A2-C8 flares 


1.88±0.21 


Feldman ct al. 1997 


Peak flux 


3.5-5.5 keV & UV 


microflares 


2.18-2.23 


Porter et al. 1995 


Peak flux 


7.7-12.5 keV 


normal flares 


1.84 a 


Hudson ct al. 1969 a 


Peak flux 


10-300 kcV 


normal flares 


1.8 


Datlowe ct al. 1974 


Peak flux 


13-600 keV 


HXR microflares 


~ 2.0 


Lin et al. 1984 


Peak count rate 


25-500 kcV 


normal flares 


1.8 


Dennis 1985 


Peak count rate 


26-500 kcV 


normal flares 


1.54 


Aschwanden & Dennis 1992 a 


Peak count rate 


>25 keV 


normal flares 


1.73±0.01 


Crosby ct al. 1993 


Peak flux 


>25 kcV 


normal flares 


1.59±0.01 


Crosby ct al. 1993 


Peak count rate 


>25 kcV 


normal flares 


1.7-1.9 


Bai 1993 


Peak flux 


>30 kcV 


normal flares 


1.8-2.2 


Bromund et al. 1995 


Stellar: 










Radiated energy, M dwarfs 


0.05-2 keV 


10 30.6 _ 1Q 33.2 erg 


1.52±0.08 


Collura ct al. 1988 


Radiated energy, M dwarfs 


0.05-2 keV 


10 30.5 _ 1Q 34.0 erg 


1.7±0.1 


Pallavicini ct al. 1990 


Radiated energy, RS CVn 


EUV 


10 32.9 _ 1Q 34.6 erg 


1.6 


Osten & Brown 1999 


Radiated energy, two G dwarfs 


EUV 


1 33.5 _ 10 34.8 erg 


2.0-2.2 


Audard et al. 1999 


Radiated energy, F-M dwarfs 


EUV 


10 30.6 _ 1Q 35.0 erg 


1.8-2.3 


Audard et al. 2000 


Radiated energy, 3 M dwarfs 


EUV 


10 29.0 _ 1Q 33.7 erg 


2.2-2.7 


Kashyap et al. 2002 


Radiated energy, AD Leo 


EUV & 0.1-10 keV 


10 31 1 - 10 33 - 7 erg 


2.0-2.5 


present work 



a quotcd in Crosby ct al. 1993 



-7- 



In statistical flare studies, the identification of weak flares close to the apparently quies- 
cent emission level becomes an ill-defined problem. Explicit detection methods discriminate 
against small flares due to overlap with larger flares, confusion between the many approxi- 
mately simultaneous small flares, or short detection times above the significance level (Hud- 
son et al. 1969; Audard et al. 2000). Limited signal-to-noise ratios add to the problem. Some 
of these complications can be overcome to some extent by fully modeling the superposition 
of a statistical ensemble of flares and comparing observable quantities between models and 
observations. An important diagnostic of this kind may be the coronal differential emission 
measure distribution (DEM) since it is determined by the ensemble of plasma packets heated 
to different temperatures. A time-integrated observation of a time-dependent heating mech- 
anism may therefore reveal a characteristic DEM that contains diagnostics for the stochastic 
flare distribution (Giidel 1997; Giidel et al. 1997; the modeling described in the latter pa- 
pers concentrated on small a with important effects due to large and long-lasting flares). 
We will derive an analytical expression for a DEM of a flare-heated corona in this paper. 
The present work and its companion paper (Kashyap et al. 2002) are, on the other hand, 
predominantly concerned with statistical light curve analysis, comparing model light curves 
with long monitoring observations obtained in the EUV and X-ray ranges. While Kashyap 
et al. (2002) focus on a detailed instrumental modeling of the EUVE DS instrument and a 
statistical analysis of the photon arrival time distribution for a sample of stars, the present 
paper concentrates on a dedicated, long observing campaign of AD Leonis using the EUVE 
and BeppoSAX satellites; two complementary methods, one presented in Kashyap et al. 
(2002), are applied, and the results are used to compute the emission measure distribution 
of a continually flaring stellar corona. 

The paper is organized as follows: §3 discusses our data reduction techniques, and §4 the 
analysis methods. In §5 we present our results, and we discuss various features in §6. Finally, 
§7 presents our conclusions. Although we follow our hypothesis that all of the observed X-ray 
emission is due to a statistical ensemble of flares, it is convenient to use the term "quiescent" 
emission for the emission level at which individual flares can no longer be separated, i.e., for 
the quasi-steady emission. 



3. DATA SELECTION AND REDUCTION 

Our target for the present investigation is the nearby dMe star AD Leo. AD Leo is 
a well-studied flare star with a high flare rate. Its spectral class is M3 V, and its X-ray 
luminosity amounts to logL x = 28.95 (L x in erg s _1 ). We use a distance of d = 4.90 pc 
(see Audard et al. 2000 for further details on the target, and for references). Its quiescent 



- 8- 



count rates in the EUVE DS, the BeppoSAX LECS, and the BeppoSAX MECS used here 
are, respectively, ~0.15, 0.105, and 0.025 cts s" 1 (these values refer to the mean of all flux 
bins in which no significant flares were evident.) 

The Extreme Ultraviolet Explorer (EUVE, e.g. Malina & Bowyer 1991) data presented 
here were obtained in several segments between 1999 April 2 and 1999 May 15 (see Table 2). 
DS Remapped Archive QPOE files were rebuilt using the euvl .8 package within IRAF. For 
our method 1 (see below) in which we used binned data, we applied primbsching (telemetry 
saturation) and deadtime corrections to all data sets, and excluded time intervals for which 
the combined correction factors exceeded 30%. To avoid fluctuations that may still occur 
within one satellite orbit window owing to residual inaccuracies in these corrections (e.g., 
due to the South Atlantic Anomaly), we decided to generally bin data to one bin per orbit 
(5663 s), which turns out to be sufficient to recognize numerous EUV flares, to resolve them in 
time, and to provide a good signal-to-noise ratio per bin. The light curve is shown in Figure 1. 
The last time interval (see Table 2) suffered from too high radiation, with correction factors 
higher than typically acceptable. This segment was not used in the subsequent analysis (for 
any of the methods). The final data set used for method 1 contained 470 bins. 

The BeppoSAX satellite (Boella et al. 1997a) pointed its Low and Medium Energy 
Concentrator Systems (LECS and MECS, respectively; Parmar et al. 1997; Boella et al. 
1997b) and its Phoswich Detector System (PDS; Frontera et al. 1997) towards AD Leo three 
times (Table 2), spanning a total time of 15 days for 270 ks of exposure time. Most of the 
observations were performed simultaneously with the EUVE observations. No significant 
signal was detected in the PDS instrument, even during flares. The three different pointings 
were similar, allowing us to merge the individual data sets into single LECS and MECS data 
sets. The cleaned and linearized event files from LECS and MECS23 (MECS2 and MECS3 



Table 2. OBSERVING LOG (1999 APRIL/MAY) 



Segment 


Instrument 


from UT 


to UT 


HJD range 






MM/DD hh:mm 


MM/DD hh:mm 


-2440000.5 


I 


EUVE 


04/02 16:09 


04/04 04:54 


11270.673-11272.204 


II 


EUVE 


04/05 00:46 


04/14 16:41 


11273.032-11282.695 


III 


EUVE 


04/17 03:32 


04/24 09:27 


11285.147-11292.394 


IV 


EUVE 


04/25 16:29 


05/04 11:41 


11293.687-11302.487 


V 


EUVE 


05/06 16:41 


05/16 05:11 


11304.695-11314.216 


I 


BeppoSAX 


05/01 06:34 


05/03 04:48 


11299.274-11301.200 


II 


BeppoSAX 


05/08 05:47 


05/10 13:10 


11306.241-11308.549 


III 


BeppoSAX 


05/12 08:39 


05/15 15:05 


11310.360-11313.628 
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2.5 



I 



II 



III 



IV 



V 



2.0 



EUVE DS 



1.5 




0.5 





O.OL 
11270 



11280 



11290 
time (HJD - 2440000.5) 



11300 



11310 



Fig. 1. — EUVE DS light curve of AD Leo (before primbsching and dead time corrections), 
obtained between April 2, 1999, and May 16, 1999. Segment V suffers from "dead spot" re- 
duction in effective area and from high radiation. The la error bars are typically ±0.01 ct s _1 
and have been plotted. 

combined file) from the SAX Science Data Center pipeline processing were filtered with the 
good time intervals that exclude the events occurring while there was no attitude solution. 
Because of the overall faster variability seen in the X-ray range and in order to maximize 
the number of available points while still retaining a sufficient signal-to-noise ratio to detect 
a large number of flares, we decided to resolve each orbit visibility interval into a few bins of 
length 200 s (a test performed with 1 bin per orbit provided compatible results that were, 
given the small number of points, ill-constrained). We thus used a total of 658 bins for the 
LECS light curve and 1363 bins for the MECS light curve. The MECS has more bins since 
it observed longer (up to about 3000 s) during each orbit, and it suffered from fewer bad 
time intervals during the on-source observations. The light curves of all three segments are 
shown in Figure 2. 

Average quiescent spectra were extracted for LECS and MECS, excluding obvious large 
flares (but still including small flares because, in our working hypothesis, there is no strict 
difference between flaring and non-flaring emission). The source was extracted inside circular 
regions of radius 8.17' and 4' for LECS and MECS, respectively. Blank sky pointings were 
used to model the instrumental and sky X-ray backgrounds. We checked the local background 
with two semi-annular regions for LECS (see Parmar et al. 1999 for full details) and two 
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0.8 
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0.2 




0.0 
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EUVE DS 



H 1 1 1 1 1- 



.'.<< .,»"' 



Ill 



SAX LECS (0.1-10 keV) 



] — i — t — — i — i — h 



SAX MECS (2-10 keV) 



11300 11301 11307 11308 11311 11312 11313 
time (HJD - 2440000.5) 



Fig. 2. — BeppoSAX MECS and LECS light curves, compared with simultaneous EUVE DS 
data. Roman numerals refer to the BeppoSAX segments (see Table 2). Note that EUVE 
DS data obtained during BeppoSAX intervals II and III are unreliable due to high radiation 
background and were not used. 
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regions perpendicular to the on-board calibration sources for MECS. We found that the 
blank-sky backgrounds are suitable for our data. We used the lemat tool in SAXDAS 2.0.1 to 
create LECS response matrices, and we adopted the standard MECS responses distributed by 
the mission (September 1997). Finally, the data were grouped with a minimum of 25 counts 
per bin. Data from 0.12 to 4 keV were kept for the LECS spectrum, while data from 1.65 to 
10.5 keV were kept for the MECS23 spectrum (see Fiore, Guainazzi, & Grandi 1999). We 
fitted the combined spectra in XSPEC (Arnaud 1996) using several isothermal collisionally 
ionized equilibrium plasmas (the so-called MEKAL code) with a fixed interstellar hydrogen 
absorption column density Ah = 10 18 cm -2 (Sciortino et al. 1999). We also introduced a free 
constant factor to the LECS spectrum to account for cross-calibration discrepancies in the 
overall effective areas. We left the Fe abundance free while we kept the other abundances at 
solar photospheric values. A three-temperature model provided a reasonable fit (Table 3). 
Finally, we derived an emission measure distribution (DEM) from the combined LECS and 
MECS data in SPEX, averaging the DEMs from the polynomial and the regularization 
methods as described in Kaastra et al. (1996a) (this is further discussed in § 6.4). 

4. ANALYSIS 

We analyzed the calibrated and cleaned data using two different methods. Since the 
present data quality does not allow us to model the run of temperatures and emission mea- 
sures of individual flares (and therefore to integrate the modeled radiative losses across com- 
plete flares), we will assume that the observed count rate is proportional to the total X-ray 
flux from the star. We will discuss the reliability of this assumption in §6.3. 



Table 3. FIT TO THE QUIESCENT BeppoSAX DATA 



logJV H a 


fcTi 


kT 2 


kT 3 


EMi/10 50 


EM 2 /10 50 


EM3/IO 50 


Fe b 


F 


x 2 /dof 


(cm- 2 ) 


(keV) 


(keV) 


(keV) 


(cm- 3 ) 


(cm- 3 ) 


(cm- 3 ) 








= 18.0 


0.18 ±0.01 


0.68 ±0.01 


1.80 ±0.08 


5.4 ±0.4 


13.6 ±0.4 


4.5 ±0.2 


0.54 ±0.03 


0.69 ± 0.03 


315.48/261 



a Fixed value 

b Abundance relative to the solar photospheric value (using values from Anders & Grevesse 1989) 
c Effective-area cross-calibration factor LECS/MECS23 
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4.1. Count Rate Distributions 

We compare the count rate distributions of the light curves with simulated data sets 
composed of a statistical flare distribution. The simulations are performed as follows: We 
first calculate a statistical power-law distribution of flare energies (equation 1) specified by 
three free parameters: the total number of flares during the simulation time (flare rate 
Rf), the power-law index a, and a lower cut-off energy E . The latter is necessary for 
practical reasons but must also physically exist in the case of a > 2. The count rate 
normalization is arbitrary and will be re-normalized later. The flares, initially defined as 
delta functions, are randomly distributed in time. The length of the simulation time interval 
was chosen to contain 5490 bins, i.e., about ten times more than the binned EUVE and 
BeppoSAX/LECS light curves and about 4 times more than the MECS light curve, i.e., we 
compare the observation with about ten (resp. four) statistical realizations of a simulation 
of equal length. 

This model light curve is then convolved with a calculated exponential decay profile 
with pre-set decay time constant r. This latter form is typical of that commonly observed 
for stellar flares. The model flares start at their peaks (i.e., the rise phase of the flare is 
very short, compatible with the observations). We note that the individual flare decay times 
vary somewhat, also due to superimposed flares that may widen the apparent flare profile. 
Experiments with a profile that was half as wide for the LECS and the MECS data showed, 
however, that the results are quite insensitive to the precise profile shape. The flare decays 
are generally considerably slower in the EUV range than in soft X-rays. Based on a detailed 
analysis with fully resolved EUVE DS data sets including ours, a decay constant of 3000 s 
proves appropriate (Kashyap et al. 2002), and we adopt this value for consistency with the 
study presented by the latter authors although we will test (as they did) variable decay times 
as well. Most rapid features in the LECS and MECS required a characteristic decay time 
scale of only 300-400 s; we adopted a decay time of 360 s. 

Since during a considerable fraction of the satellite orbit the target star is not acces- 
sible, we also simulated the effects of "windowing". The DS instrument observed AD Leo 
contiguously for typically 1520 s during EUVEs 5663 s orbit, the LECS for 1400 s and 
the MECS for 3000 s during the 5780 s BeppoSAX orbit. One anticipates that correcting 
for this effect is statistically of little importance. In particular, since the flare decays are 
nearly exponential, a flare from which the initial part is cut off behaves like a flare detected 
at peak time, with the same decay time scale but with a smaller amplitude. In this case, 
flares of a given energy statistically shift to somewhat lower energies. Since the flare decay 
constants are similar or equal for all flares, their count rates are statistically suppressed by 
the same factor, resulting only in a lower normalization of the count rate distribution, which 
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is irrelevant. Our conjecture was verified with appropriate models that cut simulated flares 
arbitrarily the same way as the real observation does. 

We then sort all bins of the light curve, both for the observed and the model light curve, 
in order of increasing count rate c to obtain a cumulative count rate distribution Af(< c) 
(number of bins with a count rate up to a given c, where we normalize Af by the total 
number of bins, i.e., < Af < 1; see Fig. 3-6). The cumulative distribution is very sensitive 
to systematic deviations in the shape of the count rate distribution as produced, for example, 
by different flare energy distributions. Since the observed and the model distributions are not 
mutually normalized in their count rates, we renormalize the (cumulative) model distribution 
such that the average count rate within a range [Afi, J\f 2 ] is the same as the average count rate 
in the corresponding Af range of the observed distribution. At this stage, noise corresponding 
to the observed Poisson noise is added to the model points, and the model data points are 
sorted again in count rate. 

The low end of the count rate distribution is inherently ill-defined. It is sensitive to 
statistical fluctuations from the superposition of numerous weak flares while at the same 
time the relative errors are largest in that range. Further, tails from a few long-decay flares 
may modulate this count rate level in time. Also, intrinsic non-flaring variability such as 
due to emerging magnetic loops or rotational modulation can introduce slow variations. We 
therefore keep the final normalization interval above zero, namely at [Afi,Af?\ = [0.1,0.5], 
and do not further consider the portion at Af < 0.2 (< 0.1 for LECS/MECS). 

For any selected a, we generated a large family of models that differ only in their flare 
rates. We selected the example that minimizes the largest vertical distance between the 
model and the observed cumulative distributions (i.e., minimal max(|A/" bs —Af mo dci\)) as our 
best fit. The traditional Kolmogorov-Smirnov test cannot be applied to the two distribu- 
tions since the normalization (one of the required fit parameters) has been derived from 
the properties of the distributions themselves. We rather simulate a sample of observations 
with defined properties (a, Rf, and E ) and a length equal to the real observation. We 
then analyze each simulated observation precisely the same way as the real observation in 
order to find statistical approximations to the confidence limits of our results. If the model 
distribution is too shallow (a too small) we find an excess of large count rates, i.e., the model 
cumulative distribution lies below the observed cumulative distribution toward larger count 
rates (see Fig. 3-6 below). Conversely, if the model distribution is too steep (a too large), 
then the model distribution lies above the observed distribution. 
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4.2. Analysis of Photon Arrival Time Differences 

Kashyap et al. (2002) present a detailed description of the second method that is founded 
on a detailed modeling of the EUVE DS detector. We provide only a brief summary. The 
basis of the method is that, as count rates rise and fall with flaring activity, the intervals 
between photon arrival times decrease and increase according to the Poisson distribution 
appropriate for the count rate at any given moment. These changes in the arrival time 
differences (St) cast a signature on the photon event list that changes according to the nature 
of the underlying source variability. Thus, for a given observation, the observed distribution 
f(St) summarizes the character of the variability of the source during that observation. A 
fixed flare distribution (equation 1) gives rise to a definite f(St) provided that the observation 
is of sufficient duration to contain a representative range of intensities, and different flare 
distributions will give rise to different arrival-time difference distributions. Note that f(St) 
is not sensitive to the actual temporal locations of the flares, but rather depends only on the 
stochastic ensemble described by equation 1. The power-law index, a, that best describes 
the observed light curve is determined by comparing the observed arrival time difference 
distribution f b s (St), with simulated distributions f S i m (Si). 

This method was developed to deal easily and rigorously with the windowing inherent 
in low Earth orbit observations (with a typical observing time of 30 minutes for each ~ 
90 minute orbit), such as those obtained by EUVE, and to allow for proper treatment 
of telemetry saturation ( "primbsching" ) and deadtime effects that can introduce somewhat 
variable corrections during one orbit. Since we deal here with photon lists, a proper treatment 
of the corrections continuously in time is important. The advantage of the method is that 
it operates directly on the observed photon event list, avoiding the need for time binning. 
Simulated event lists can be windowed in exactly the same way as the observed event list, 
and primbsching effects can be applied to the simulated events, censoring them in the same 
stochastic way as the observed events. 

We assume here that the observed light curves can be described by the sum of a flaring 
component, with a power-law frequency distribution of flare energies as described by equa- 
tion 1, and a constant component. The flaring and constant components are described by the 
flare rate Cf, and a "quiescent background" rate Cb, respectively. As for the first method, 
the flares are assumed to be impulsive events whose count rates decay exponentially with a 
time scale of a few thousand seconds. The exact decay time can be varied to best match any 
observed flares or the observed and synthetic photon event lists themselves. In practice, we 
have found that results obtained for decay times in the range 1000—5000 s are typically very 
similar because the extensive overlapping of flare events tends to decrease the sensitivity of 
the results to this parameter (Kashyap et al. 2002); the final results presented here assume 
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a decay time of 3000 s, identical to the value found and applied by Kashyap et al. (2002). 
For specified values of Cf and Cb, a synthetic light curve corresponding to the entire interval 
covered by the observed light curve can then be realized through a Monte Carlo algorithm, 
assuming a random distribution of flares in time but subject to the power-law frequency dis- 
tribution of total energies. This light curve is then windowed by the observed photon event 
"good time intervals" and a synthetic event list is derived through a Poisson realization of 
the resulting light curve. The synthetic event list is then pruned by discarding photons at 
a rate corresponding to the observed Primbsch factor. The remaining set is identical in its 
"instrumental characteristics" to the observed data and the observed and synthetic event 
lists can be compared directly. 

Observed photon arrival time differences are compared to those synthesized across a 
grid of the parameters a, Cf and Cb using the x 2 statistic to compare f S i m (5t) with f b s (St). 
A number of simulations (typically ~ 10) are carried out and the median value among the 
resulting x 2 are used to compute the likelihood of obtaining the observed data for the given 
set of parameters {a,Cf,Cb}- In the case of AD Leo, we have examined each of the three 
EUVE observation segments II, III, and IV independently (Table 2), treating them as three 
different observations, as well as treating the whole sequence II-IV at once. This enabled us 
to examine the degree of consistency of our derivation of a from segment to segment. 



5. RESULTS 

5.1. Parameters and Tests 

Our simulated light curves were generated assuming an exponential flare decay time 
r independent of the flare energy or amplitude. If the decay time increases systematically 
with the total radiated energy, e.g., r oc E 13 , then the flare amplitudes increase less than 
proportionally with the energies, which may be interpreted as a steeper energy distribution. 
However, increasing the decay time of larger flares produces more bins at large count rates 
(but only an equal number of additional bins at low count rates) , and this effect counteracts 
the apparent steepening of the energy distribution. 

Observationally, a dependence of the decay time on energy is marginal. Aschwanden et 
al. (2000) investigated scaling laws from solar nanoflares to larger flares, covering 9 orders 
of magnitude in energy. While most of the geometric and physical parameters exhibit a 
strong scaling with the flare size, Aschwanden et al. report that the time scale (radiative or 
conductive) does not depend on the flare size. In fact, as stronger flares tend to show larger 
electron densities and higher coronal peak temperatures (§3.1 in Aschwanden et al. 2000), 
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they should radiatively or conductively decay faster than small flares, which is not observed. 
The range of decay times (their Table 2) is much larger than any possible trend between 
nanoflares and large flares. Feldman et al. (1997) measured soft X-ray FWHM durations of 
a large sample of flares over three orders of magnitude and found no trend. The differences 
between flare decays are due to scatter within classes of flares of about equal energy. Shimizu 
(1995) investigated energies and durations of small solar active-region transients and again 
found no clear trend in the durations over four orders of magnitude in peak count rate. We 
inspected the best defined flares in our data although superpositions of flare profiles may 
introduce ambiguity. The MECS data are best suited for such an analysis since they provide 
the longest intervals of visibility per orbit, a high flare-to-quiescent emission contrast, and 
a time bin resolution of 200 s. We found no trend for a correlation between flare amplitude 
and decay time; the scatter in the decay time itself dominates. The same holds true for 
the other data sets. However, the sample of selected large X-ray flares of Pallavicini et al. 
(1990) from different stars shows a weak trend if four orders of magnitude in radiated energy 
are included, but the total variation is no larger than the scatter in duration at a given 
energy. For individual stars, there is no clear trend. Their data sample can be best fitted 
with a relation r oc E - 25 . Given that our dynamic range (ratio between largest to smallest 
explicitly detected flare count rate) is of the order of 10—20 for binned data, this effect may 
statistically influence the results. We therefore tested the DS, LECS, and MECS analysis 
by introducing variable decay times as specified above. 

We have used two realizations of our analysis for the binned EUVE data (method 1) 
that show a well-developed quiescent emission. First, we applied it to the data from which 
the (constant) quiescent level was subtracted. In this case, we test whether we can find a 
model distribution that is compatible with the observed flare emission. Second, we applied 
the method to the complete data including the quiescent emission. Since the latter, in our 
model, should be the superposition of unresolved small flares, the addition of a quiescent level 
simply corresponds to the extrapolation of the power law to lower energies, and we expect 
that the power-law index does not significantly change. This second realization also allows 
us to derive the lower cut-off energy Eq required to explain the quiescent emission. The 
comparison between the two realizations could potentially reveal a basic difference between 
quiescent and flaring emission. If the flare energy distribution does not steadily continue 
toward more numerous small flares that eventually merge with the quiescent level, then the 
first realization would be subject to a cut-off energy possibly above the quiescent level (bi- 
modal flux distribution), and the derived values for a may differ. These two extreme cases 
will further be discussed below. The contrast between flares and the quiescent emission is 
much larger in the BeppoSAX data; we treated only the complete data sets in these cases for 
the primary results but also subtracted the quiescent level for an assessment of the confidence 
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ranges (see below). 

5.2. Basic Findings 

We first discuss results obtained using method 1 under the assumption of constant 
exponential decay times r for all flares. Figure 3 shows results for the EUVE light curve of 
AD Leo (segments I— IV) from our first method. The top figure represents the light curve 
binned to one point per orbit (5663 s). Here, a constant quiescent count rate of 0.13 ct s^ 1 
has been subtracted, corresponding to the lower envelope of the light curve. The x-axis 
("time") gives the sequential bin number, with observing gaps at 30, 200, and 340 fully 
considered. The alternative analysis in which the initial large flare was excluded used only 
bins above bin no. 70. 

The three lower rows then illustrate, in this order, the best fit found in the analysis (i.e., 
optimum a and optimum flare rate), a selected case with too low a (but again optimum flare 
rate), and similarly a case with too high a. In each case, the figure on the left shows an 
extract of the simulated and normalized light curve with noise added, while the figure on 
the right shows the respective cumulative distributions. In the latter, the solid line shows 
the observed cumulative count rate distribution Af(< c) while the dashed line illustrates the 
normalized count rate distribution from the model. The two dotted vertical lines indicate 
the range used to normalize the model distribution to the observation. Our statistical fit 
criterion (minimizing the largest vertical distance between the two distributions) was applied 
only above the dashed horizontal line (above M = 0.2 or 0.1). The dotted function close to 
and around the horizontal M = line illustrates the difference "model - observation" for the 
cumulative distributions. Evidently, the optimum index is a ~ 2.2 in this case, while values 
as low as 1.9 or as high as 2.7 are considerably worse. 

Figure 4 shows the equivalent analysis but without subtraction of a quiescent level. The 
optimum value is formally found at a — 2.1 although the cumulative count rate histogram 
suggests a somewhat higher a. The fit is generally somewhat poor, but a test in which the 
quiescent level was reduced by 50% still produced a = 2.1. Lastly, Figures 5 and 6 show 
the results from the analysis of the LECS and MECS light curves, respectively. Only the 
observed light curve, an extract of the optimum model, and the cumulative distributions are 
shown. We find optimum values of a = 2.4 and 2.2, respectively. 

The statistical confidence ranges derived from our sample of simulated observations are 
as follows. Approximately 90% of best-fit solutions for the EUVE DS simulations (performed 
for a = 2.0 and a = 2.2) were found symmetrically within a range of ±0.1 centered at the 
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Fig. 3. — Examples of statistical flare simulations from method 1 for different power-law distribu- 
tions. The EUVE DS data set of AD Leo is shown. Only emission exceeding quiescent level has 
been modeled. The early large flare is included. Three simulated best-fit examples for different a 
are shown. Top figure: Observed light curve. Second line: Optimum case; a = 2.2. Third line: 
Too shallow distribution with a = f.9. Bottom line: Too steep distribution with a = 2.7. For 
2nd, 3rd, and bottom figure panel: Left: Simulated light curve (extract, normalized). Dashed line 
marks the lower limit to count rates that were considered for optimizing the fit. Right: Cumulative 
count distribution for data (solid) and model (dashed), and difference (dotted, around zero line). 
The maximum difference (vertical bar between observation and model) is minimal for the best 
fit. The two vertical dotted lines mark the interval for the model count rate normalization. The 
horizontal dashed line marks the count rate level above which the differences between observation 
and model were considered. 
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Fig. 4. — Similar to Fig. 3 (method 1), but all emission, including quiescent emission, 
has been modeled. Top panel: Observed light curve. Second panel: Optimum case; 
a — 2.1. Third panel: Too shallow distribution with a = 2.0. Bottom panel: Too steep 
distribution with a — 2.7. 
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Fig. 5. — Similar to Fig. 3 (method 1), for LECS. Only optimum model is shown, with 
a = 2.4. Gaps between observation segments and between orbits are not retained. 
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Fig. 6. — Similar to Fig. 3 (method 1), for MECS. Only optimum model is shown, with 
a = 2.2. Gaps between observation segments and between orbits are not retained. 
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correct a. We thus adopt Ao = ±0.1 as our 90% confidence limits for the EUVE DS 
results. The lower signal-to-noise ratio of the BeppoSAX data constrained the simulation 
results less well, probably due to the large dominance of the noise fluctuations at low count 
rates. We therefore repeated the analysis (both for the real and the simulated observations) 
constructing cumulative distributions from those data only that exceed the quiescent level 
by lcr, thus selecting predominantly the visible flares only. For the LECS, a = 2.4 was 
confirmed, while for the MECS, we found optimum fits for a = 2.0. From the sample of 
simulated observations, we again find a 90% width of about ±0.1, although a is slightly 
underestimated by 0.1 on average. This effect is probably due to selectively choosing the 
larger count rate bins. We conservatively adopt a combined (statistical and systematic) 
uncertainty of ±0.2 for the LECS and MECS (90%). 

The following systematics are evident (Table 4): (i) In all cases, we find acceptable 
values for a at or above 2, with a — (2.1 — 2.3) ±0.1 for the DS data, a = 2.4 ± 0.2 for 
LECS and a = (2.0 - 2.2) ± 0.2 for MECS. When we excluded the large flare in the DS light 
curve, the distribution steepened by Aa ~ 0.1 (Table 4), an effect that is to be understood 
as follows: The presence of a large flare will tend to give a smaller a because comparatively 
less power will be contained in smaller flare events. Conversely, the selective elimination 
of the population of the strongest flares steepens the distribution although a genuine single 
power law cannot be retained if too large a fraction of the distribution is eliminated. The 
long decay of the large flare clearly also biases the count rate distribution toward a lower 
a (adding a large number of bins at high count rate levels). Overall, thus, the DS and the 
LECS ranges for a agree. The MECS range appears to be systematically lower. We will 
discuss this effect further in Section 6.3. 

The results from the second method, based on photon arrival statistics, are illustrated in 
Figures 7a— 7d, where the derived probabilities of power law indices a matching the observed 
index are plotted as a function of a. The results are also summarized in Table 5. The first 
three figures illustrate the derived probabilities for the three segments II, III, and IV treated 
separately, while the fourth shows the results of the analysis of the whole data set II-IV 
treated as a single observation. The important feature of all these figures is that the most 
probable value of a is again always greater than 2 but less than 2.3, based on 90% or 95% 
confidence intervals. Segment II appears to have an optimum index that is slightly lower 
than that of the last two segments. This is caused by the large flare, which contains a 
significant fraction of the total observed counts in that segment (see above). The confidence 
intervals for the last two segments are remarkably similar, indicating 2.1 < a < 2.3. The 
most probable value of a based on all three segments (Figure 7d) is a = 2.2. The values 
found here are thus in excellent agreement with results from the first method. 
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Table 4. a VALUES FOR AD LEO FROM METHOD 1; MINIMUM MODEL FLARE 

ENERGIES AND LUMINOSITIES a 



Data set, 


bin 


a 


Minimum flare 


Minimum flare 


Minimum peak 


Detection limit 


Detection limit 


segment 


size 


(90% error) 


counts 


energy 


luminosity 


flare energy 


peak luminosity 




(s) 






(erg) 


(erg s- 1 ) 


(erg) 


(erg s" 1 ) 


DS I-IV b 


5663 


2.2(0.1) 


62 


1.7 x 10 31 


5.8 X 10 27 


3 x 10 31 


1 x 10 2S 


DS I-IV C 


5663 


2.3(0.1) 


60 


1.7 x 10 31 


5.5 x 10 27 


3 x 10 31 


1 x 10 2S 


DS I-IV d 


5663 


2.1(0.1) 


0.5 


1.4 x 10 29 


4.8 x 10 25 


3 x 10 31 


1 x 10 28 


DS I-IV C 


5663 


2.3(0.1) 


3.5 


9.6 x 10 29 


3.2 x 10 26 


3 x 10 31 


1 x 10 2S 


LECS I III 


200 


2.4(0.2) 


0.8 


3.8 x 10 29 


1.0 x 10 27 


1.4 x 10 31 


4 x 10 2S 


MECS I III 


200 


2.0-2.2(0.2) 


0.4 


1.1 x 10 30 


2.9 x 10 27 


3.6 x 10 31 


1 x 10 29 



a For given time bin size and for the optimum case; minimum flare counts/energy /peak luminosity refer to smallest flares used in the 
simulation with optimum KS test result. Detection limit refers to 3<r levels for the respective bin size. 

b quiescent level subtracted, all data 

c quiescent level subtracted, large flare excluded 

d quicscent level included, all data included 

e quiescent level included, large flare excluded 



Table 5. a VALUES FOR AD LEO FROM METHOD 2 (DS) 



Segment 


Most probable 


95% Confidence 


DS 


a 


Interval 


II 


2.07 


2.00-2.13 


III 


2.22 


2.11-2.31 


IV 


2.25 


2.13-2.30 


II-IV 


2.19 


2.14-2.23 
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AD Leo: 99apr05-99opr09 (cf=2.07^^ ) 



0.10 - 



0.01 - 




1.90 2.00 2.10 2.20 2.30 2.40 

AD Leo: 99opr1 7-99opr24 (a = 2.22*^J ) 



0.100 r 



0.010 r 



0.001 - 




1.90 2.00 2.10 2.20 2.30 2.40 

AD Leo: 99opr25-99moy04 (q 2.2b>^ ) 




2.40 



0.01 r 



Fig. 7. — The derived probability of the power-law index a matching that of the EUVE DS 
observations of AD Leo based on the distributions of photon arrival time differences (method 
2). The most probable value of a corresponds to the peak value of each distribution and is 
printed together with the corresponding 95% confidence intervals at the top of each figure. 
From top to bottom: a) Segment II, including the large flare, b) Segment III. c) Segment 
IV. d) Segments II-IV combined. 
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We further investigated the case for an exponential decay time r oc E 13 = E - 25 . In the 
first method, sufficiently small flares will have decay times much smaller than the bin length 
and thus the method is sensitive only to their total energy contribution, but not to their 
decay time. Given our relatively coarse binning, we thus expect to find rather stable results 
for a. This is confirmed, with deviations of no more than ±0.1 in a. The second method is 
much more sensitive to details in the light curve and could, in principle, find both power-law 
indices (5 and a. In practise, however, there is an acceptable family of solutions (a,/3). The 
acceptable a increases with increasing (3, reaching values up to 2.6—2.8 for (3 = 0.25. We 
interpret this as being due to a larger time occupation by relatively large count rates, i.e., 
short photon arrival time differences, if the energetic flares are stretched in time relative to 
the small flares (despite the corresponding decrease in peak count rate). The model light 
curve thus appears to be "too hard" , and a higher a is required to match the statistics of the 
observation. If we require agreement with method 1, then /3 — and a = 2.2, compatible 
with the absence of a detectable trend in the decay times. Even without constraining method 
2, the lowest possible a range (for the lowest acceptable (5 — 0) is around 2.2, i.e., values 
below 2 remain excluded. 



6. DISCUSSION 

6.1. The Validity of Power Laws 

The present analysis was primarily motivated by the finding that the occurrence fre- 
quencies of solar flare energies are distributed in power laws. Although some statistical 
models support such distributions (e.g., Lu & Hamilton 1991; Vlahos et al. 1995) their phys- 
ical cause is a matter of debate. Indeed, there is little reason to assume that the same power 
laws strictly hold for all classes of flares and for all coronal regions (active regions, magnetic 
areas above the network, quiet regions, etc). A simple estimate, for example, suggests the 
presence of a high-energy cut-off: An average coronal magnetic field of strength 100 G in 
a half-spherical active region of radius 10 10 cm contains 8 x 10 32 erg of magnetic energy. 
Stronger flares (as often observed in active stars) either require stronger magnetic fields or 
a larger volume, but both are constrained by the maximum magnetic field strength avail- 
able in the photosphere and, respectively, by the volume of active regions (that are related 
to the size of bipolar magnetic regions in the photosphere and by about two coronal scale 
heights; Serio et al. 1981). A high-energy cut-off is also suggested from avalanche model 
simulations (Lu & Hamilton 1991; Lu et al. 1993) and evidence for its existence has been 
found in observations of small solar active regions (Kucera et al. 1997), and in investigations 
of time-dependent a values for solar flares (Bai 1993; Bromund, McTiernan, & Kane 1995). 
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On the low-energy side, the power-law index may change due to physical changes in the 
energy release (Hudson 1991) and this seems to be the case for solar microflares as opposed 
to the larger flares (Krucker & Benz 1998). A low-energy cut-off is also required to confine 
the total radiated power if a > 2 (see §1). 

We do not construct the flare distributions explicitly with our methods, but use model 
light curves for statistical comparison. It is, however, clear that single power laws describe 
the observed light curves acceptably well. Although flare energies detected by these ob- 
servations are comparable to medium-to-large solar coronal flares, we find a values that 
are markedly larger than those reported for equivalent solar flares, but that resemble those 
recently reported for solar microflares with ^ 10 6 times smaller energies. Typical activity 
indicators for active stars such as ours (like L x /L ho \, or the surface X-ray flux) are of order 
1000 times larger than the Sun's. It thus appears that the regime of steep power-law indices 
(a > 2) is shifted upwards in energy on active stars, and that we see an equivalent popula- 
tion of flares at larger energies. We find power-law indices around 2.0-2.5 and our methods 
exclude values below a = 2 for the softer DS and LECS light curves. Furthermore, the 
smallest flare energies that we require to fully model the observations lie considerably below 
the flare detection limit in the binned light curves (Table 4) although our second method 
may reach a sensitivity close to such values (Kashyap et al. 2002). Like the solar microflares, 
our population of flares may thus be sufficient to energize the complete corona, including 
the quiescent emission, if the power-law is extrapolated to flares with radiated energies of a 
few times 10 29 erg (Table 4). Such energies correspond to relatively moderate flares in the 
solar context. These conclusions fully support the findings by Audard et al. (1999, 2000). 



6.2. The Minimum Flare Energies and Quiescent Emission 

A lower cut-off to the flare energies is required if a > 2. The cut-off does not imply that 
lower-energetic flares do not exist. But it implies that their occurrence rate cannot follow the 
extrapolation of the power-law found at higher energies but must be considerably smaller, 
effectively introducing a cut-off below which flares contribute little. Our simulated model 
light curves were calculated assuming such an energy cut-off. We determined the cut-off 
energy after renormalization and thus found the minimum number of counts of any of the 
simulated flares. The values are reported in Table 4, column 4. To convert the total number 
of counts to total (radiated) flare energy, we computed the count-to-energy conversion factor 
as follows: A 3-temperature model (reported in Table 3) determined from the LECS and 
MECS data by spectral fitting of the quiescent emission was convolved with the response 
matrices of each of the three instruments. We thus obtained the expected total count rates 
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in the detectors, and by integrating the complete spectrum from 0.01 keV to 50 keV we 
computed the total coronal luminosity. We thus find that one detected (quiescent) count 
corresponds to 2.8 x 10 29 erg for the EUVE DS, to 4.8 x 10 29 erg for the BeppoSAX LECS, 
and to 2.6 x 10 30 erg for the MECS. From this, we obtain the minimum energy of flares used 
for the simulation, as shown in column 5 of Table 4. In other words, extrapolating the power- 
law distribution of flares down to the reported energies is necessary and sufficient to explain 
the complete observed quiescent flux. The LECS provides a lower limit of 4 x 10 29 erg. For 
the derived flare time profile, this energy corresponds to a peak luminosity of approximately 
1 x 10 27 erg s _1 in the combined X-ray and EUV ranges. Such flares correspond to small 
solar flares. Similar values hold for the MECS and DS data (Table 4, columns 5 and 6). 

We note that the actually detected flares in the binned light curves (method 1) typically 
exceed these levels by ~2 orders of magnitude. A small flare reaching a peak count rate 3a 
above its pre-flare level corresponds to 3 x 10 31 erg in the DS, 1.4 x 10 31 erg in the LECS, 
and 3.6 x 10 31 erg in the MECS (Table 4, column 7). The apparently quiescent level is thus 
composed of flares between the cut-off limit and the 3a count rate detection limits. In terms 
of peak luminosity, this interval covers the ~ 5 x 10 25 — 10 29 erg s -1 range (see Table 4, 
columns 6 and 8, for details). Assuming an average radiative loss function (radiative energy 
loss per unit EM) of A = 2 x 10~ 23 erg cm 3 s _1 (appropriate for T = 5 — 40 MK) we find that 
the peak EMs of these flares are approximately 2.4 x 10 48 — 5 x 10 50 cm -3 for the DS quiescent 
level, 5 x 10 49 - 2 x 10 51 cm" 3 for the LECS quiescent level, and 1.5 x 10 50 - 5 x 10 51 cm" 3 
for the MECS quiescent level. From Figure 2 of Feldman et al. (1995) we estimate that the 
peak temperatures of these flares are 15-30 MK for the DS, 25-33 MK for the LECS, and 
27-37 MK for the MECS (after Aschwanden 1999, these temperatures are smaller by a factor 
of ~ 1.5). If we consider that most of the flares in the power-law distribution are close to 
the lower end of the temperature intervals, we see that the LECS quiescent level is primarily 
composed of flares that reach no more than about 20 MK at peak but mostly reside at lower 
temperatures. The MECS effective area shows a steep gradient below 25—30 MK, thus most 
of the LECS quiescent emission is suppressed in the MECS. The 3-T fit explicitly shows that 
the bulk quiescent plasma is at temperatures around 10 MK (Table 3). 

The two BeppoSAX light curves indeed look qualitatively different (Fig. 2). Although 
the same flares are present in both light curves, the LECS shows appreciable quiescent emis- 
sion while the MECS does not. The contrast between strong flares and low- level emission is 
much stronger in the latter. This difference is a consequence of different temperature sensi- 
tivities of the detectors, the MECS being insensitive to low-level emission from presumably 
cooler plasma. To quantify this hardening effect, the count rates were first normalized to the 
average quiescent count rate level. Fig. 8 (left) illustrates the relation between LECS and 
simultaneous MECS count rates for 658 bins of 200 s each (see Fig. 2) along the complete 
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BeppoSAX light curve. It shows a pronounced deviation from proportionality: the harder 
MECS emission increases faster than the corresponding softer LECS emission, i.e., the emis- 
sion hardens with increasing overall count rate. The best-fit power-law in this figure has an 
index of 1.45 ± 0.02. 

We thus find a natural explanation for the largely differing contrast between flares 
and quiescent emission in the BeppoSAX detectors (Figure 8 left) in the framework of the 
stochastic-flare heating hypothesis. We next need to investigate whether the larger flares 
themselves show a trend for hardening with increasing count rate. We emphasize again that 
we measure hardening by using two different effective area curves for two detectors, not by 
applying any spectral analysis. 

There are two possible contributions that induce the non-proportionality between nor- 
malized LECS and MECS count rates in Fig. 8 (left): i) There is an intrinsic hardening for 
the larger flares, as suggested by the Feldman et al. and Aschwanden et al. relations. Thus, 
if the detector effective area curves vary differently across the relevant temperature range for 
LECS and MECS, we expect to see hardening signatures, i.e., non-proportionalities in Fig. 8. 
ii) The two detectors have, as argued above, different relative sensitivities to the quiescent 
emission compared to the (presumably hotter) large flares. To compare the flare count rates, 
the quiescent levels should be subtracted. We thus subtract a baseline quiescent count rate 
(0.105 ± 0.010 cts s- 1 for the LECS and 0.025 ± 0.005 cts s" 1 for the MECS). Fig. 8 (right) 
shows the relation between the residual count rates for LECS and MECS after the optimum 
subtraction. The slope of the best-fit line is 1.10loo5> i- e -> there is still some hardening effect 
in the flare count rates. We show in the next subsection how the hardening of larger flares 
can explain the shallower flare energy distribution found for the MECS (compared to the 
LECS results). 



6.3. Dependence on Spectral Range 

The range of best-fit a agrees for all tests performed on the data from the EUVE DS 
and from the BeppoSAX LECS (2.1 < a < 2.4). There seems to be a systematic shift of the 
MECS results relative to those from LECS, by approximately Aa = 0.2 — 0.4 to lower a, 
i.e., the MECS distribution is shallower although this effect is only marginally significant. 

For an isothermal plasma, a harder energy spectrum implies a higher plasma tempera- 
ture. Feldman et al. (1995) reported a correlation for solar (plus a few stellar) flares between 
flare peak temperature T and flare peak emission measure EM (which scales approximately 
with the flare peak luminosity and, for constant decay times, with the total radiated flare 
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Fig. 8. — Correlation between simultaneous LECS and MECS count rates measured in 658 
bins of 200 s duration each, for the complete light curve. The count rates were normalized 
with the quiescent values, and only normalized values exceeding unity are plotted. Left: All 
data included. Best-fit slope of the regression fit (dashed, using weights calculated from error 
bars): 1.45±0.02. - Right: Quiescent count rate level was subtracted before normalization. 
Best-fit slope of the regression fit (dashed): l.lOloios- 
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energy). The relation is nearly exponential but can be reasonably approximated also by a 
power law over a limited range of temperatures, 

EM = aT fc [cm" 3 ] (3) 

with a si 2 x 10 13 cm~ 3 K~ 6 , b m 5 ± 1 in the range of T = 5 — 30 MK, and T measured in 
K. Aschwanden (1999) reports a power-law dependence between T = 1 — 20 MK with b m 7. 
The X-ray luminosity in general can be expressed as 

L « EM A(T) = / EM T" (4) 

with (j) ~ 0.3 over the above temperature range (for broad-band X-ray losses as derived in 
XSPEC or SPEX). 

The count rate measured in a detector depends on the effective area curve and the inci- 
dent X-ray spectrum which in turn depends on the plasma temperature T and the emission 
measure. We approximate the dominant X-ray emission at any one time as being emitted by 
an isothermal plasma of temperature T, with equation 3 satisfied, i.e., we identify any count 
rate at any given time with an individual flare at its peak and thus use its peak EM = EM 
and peak T = T for this estimate. Because we cannot perform time-resolved temperature 
analysis with the present data, we approximate the true flare energy distribution with the 
observed count distributions (i.e., c oc E; the indices £ and m stand for LECS and MECS, 
respectively): 



dN 

dc£ 
dN 



(5) 

oc c~ a - (6) 



where the right-hand sides hold for the hypothesis that the energy distributions (i.e., the 
count rate distributions) are power laws. This approximation is acceptably good for out- 
standing flares, i.e., count rates well above the quiescent level, but breaks down near the 
quiescent level where numerous small flares may overlap. 

We describe the dependence of the LECS and MECS count rates q and c m on temper- 
ature T and luminosity L with power-law approximations 

d = ce, LTT (7) 
c m = c rnfi LT^ (8) 

where c^o and c mi o are detector-related constants and 7^ and 7 m describe the temperature 
sensitivity of the detector count rate. The 7 values must be determined from the detector 
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effective areas. The expression £p = cz/L = c^T 11 (similar for m) gives the observed count 
rate per unit luminosity, i.e., the efficiency of the detector to record a given luminosity in 
terms of a count rate, as a function of temperature. 

We have folded a number of theoretical isothermal model spectra (in collisional equi- 
librium, as modeled in XSPEC [Arnaud 1996], and SPEX [Kaastra et al. 1996a]) with the 
detector response matrices for the DS, LECS, and MECS detectors to describe £^ m (T). The 
count rates were derived across the range of sensitivity of the detector, and the L values have 
been evaluated for each T by integration of the model spectrum from 0.01—50 keV assuming 
a unit emission measure. The results are shown in Fig. 9. The efficiencies of the LECS and 
DS detectors are rather flat (7^ ~ 0) at least in the region of interest, i.e., within approxi- 
mately 0.5 — 5 keV. The LECS is least sensitive to changes in T and therefore best recognizes 
flares of different temperatures with the least bias. On the other hand, the efficiency of the 
MECS is a strong function of T below ~ 3 keV, with a rapid drop-off toward emission from 
cooler plasmas. Because of the dependence in equation 8, this implies that weak flares are 
suppressed in the MECS, and strong flares are enhanced. This is equivalent to an apparent 
decrease of a compared to the initial energy distribution. 

To estimate this effect, we set 7^ = (for the temperature range of interest). The 
number distribution of flares (i.e., counts) in LECS count rate, dN/dce, is therefore an 
unbiased approximation to the true flare radiative energy distribution dN/dE since E oc L. 
Since 

' '1,0 j 7m _ C m,0 / \ _ C m ,Q / Q j ^ 



Q Q,o Q,o \ a f J Q,o 

(using equations 3 and 4) we obtain 

c m occJ m/( ^ )+1 . (10) 

Therefore, 

dN = dl^dcrn ^ dN_ crtm /{b-<t>) ^ 
dee dc m dc e dc m e 
With equations 5 and 6 we find that 

(1m . -1 \ 7m /-, n \ 

According to § 6.2, two components may be responsible for the apparent decrease of a; we 
consider them to be limiting cases, (i) The hardening in Fig. 8 (left) could be attributed 
to the intrinsic spectral hardening toward larger energies of the visible flares only. We 
found c m w 4- 45±0 - 02 , hence 7 m /(6 - <p) + 1 = 1.45 and therefore 7 m /(6 - <p) = 0.45, so 
that a e = lA5a m — 0.45. We find the highest confidence for ap ps 2.4, hence we expect high 
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Fig. 9. — Efficiency (count rate/luminosity) as a function of (isothermal) plasma temperature 
for the EUVE DS, and the BeppoSAX LECS and MECS detectors (arbitrarily normalized). 
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confidence for a m 1.97±0.02 (not including possible systematic errors of ±0.1 as discussed 
above). 

ii) Preceding subtraction of the quiescent count rate level both for LECS and MECS 
suggests (Figure 8 right) c m ps c i- 10 (+ - 08 '-°- 06 ) ) h ence 7 m /(& — 0) + 1 = 1.10 and 7 m / (b — <f>) = 
0.1, so that ae = 1.10a; m — 0.10. We thus expect high confidence for a m 2.27 ± 0.08. 

The expected difference a;-a m p=j 0.13—0.43 agrees with the measured difference of 0.2 — 
0.4. We conclude that the hardening of flares with increasing total energy significantly affects 
the determination of a. While the LECS data provide, thanks to the flat efficiency curve, 
a relatively unbiased measure of the flare energy distribution (in terms of emitted energy), 
the MECS data consistently result in a harder appearance of the light curve, especially if 
the quiescent emission is included. We thus recognize the MECS result as biased by detector 
properties. As far as the DS results are concerned, they agree well with the LECS results, 
which is consistent with the relatively flat efficiency curve of the DS (Fig. 9). 

Similar temperature bias has been discussed for solar coronal flare statistics derived from 
emission line flux ratios in the EUV range covering a narrow temperature range (Aschwanden 
& Charbonneau 2002; Aschwanden & Parnell 2002). One may extend this analysis to solar 
HXR data that have traditionally been used as a diagnostic of the total flare energy release 
(Lin et al. 1984). If the production of hard X-rays is more efficient in strong flares (the "Big 
Flare Syndrome", Kahler 1982), then hard X-ray flare energy distributions are shallower 
than those constructed from soft X-ray or EUV data. There are at least indications that 
very small, low-temperature flares on the Sun are poor in radio emission, i.e., they produce 
high-energy electrons (important for the HXR production) less efficiently than larger flares 
(Krucker & Benz 2000). 

6.4. The Emission Measure Distribution 

We will now investigate effects of our hypothesis in the limiting case that all of the 
observed X-ray and EUV emission is due to superimposed stochastic-flare emission. In that 
case, we can use our a values to develop a crude model for the time-averaged EM distribution 
of the observed corona. We will estimate the amount of EM at a given temperature produced 
by a flare of any energy during its decay phase, and then integrate over the flare energy 
distribution. The time a flare spends at a given temperature is used as a weighting factor. 

From equations 3 and 4 we obtain a relation between the flare peak temperature T and 
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its peak luminosity L , 

/ T \ 1/(6-*) 

T.-(0) . (X3) 

Note that for decay times r independent of E, the peak luminosity follows the same power- 
law as the total radiated energy: 

§- = kr-^L- (14) 

(since E = L t). However, we will also investigate the case in which r varies with the flare 
energy, namely 

r = tqE 13 (15) 



(where (3 > is assumed, and r is a constant adjusted to the larger detected flares) in which 
case 

" V '' V '"-' ^'-»«'-" (16) 



dL dE dL 

where the constant k! = kr^ 1 ^ /(l — (3) > as long as (3 < 1 (which can be reasonably 
assumed). Since we neglect the short rise time of the flare, our flare light curves are described 
by their exponential decay at t > 0, 

L(t) = L e-' /T . (17) 

From hydrodynamic modeling, theory and observations, it is known that during the flare 
decays T oc n**, where n is the plasma density (Reale et al. 1993). The parameter ( is usually 
found between 0.5 and 2 (Reale et al. 1993). A value of ( = 2 holds if the heating source 
abruptly turns off at flare peak, and the flare cools freely via radiation and conduction. A low 
value of C indicates sustained heating during the decay phase. For example, for ( « 0.5 — 0.7 
the decay time of the flare emission is 2—4 times slower than predicted for a freely cooling 
magnetic loop (Reale et al. 1997). Solar observations of moderate flares with SMM show ( 
between ~0.5-2 (Sylwester et al. 1993), although Reale et al. (1997) find a predominance of 
values around ~ 0.3 — 0.7, i.e., flares with sustained heating. For larger stellar flares, there 
is much evidence for small ( as well: Reale & Micela (1998) find a flare decay that is > 2 
times slower than for a freely cooling loop; Favata et al. (2000a) find ( = 0.56 ± 0.04 for a 
flare on EV Lac, and Favata et al. (2000b) report ( = 0.48 ± 0.06 for a flare on AD Leo; 
Gudel et al. (2001b) find ( = 0.95 ± 0.15 for a moderate flare on AB Dor. Since, under the 
assumption of constant volume, the emission measure scales with n 2 , we have 

L = (^) a \ m 

T VEMq / k ' 
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From equations 4, 17 and 18, we obtain the temperature time evolution 

T(t) = Toe-'/^/C-^]) ( 19) 

where we require < 2/( for a temperature decay (this is usually fulfilled as ~ and 
2/C > 1)- The emission measure EM T at temperature T for this particular flare is, from 
equations 18 and 3, 

EM r = aT fe " 2/c T 2 ^ (20) 

and this expression must be weighted with the fractional duration (out of the total ob- 
serving time P) during which the flare resides within the logarithmic temperature interval 
[InT, InT + dlnT], i.e., dt/P. The emission measure contribution EMy thus counts in full 
if it is constantly present during the observing time. The duration dt is obtained by dif- 
ferentiating equation 19 (we only require the absolute value of rflnT), so that the weight 
is 

dt 2r (\ 



p = p{rV ldlan (21) 

Multiplying equations 20 and 21 and substituting all expressions, we find the weighted 
contribution of a flare with peak luminosity L to the emission measure at temperature T, 



5(EM T ) = 



2ar 1/(1 ^ } fl 



X 2 

(6-2/0/(6-0) 



X 



t^j L^T^\d\nT\ 
= cLf-WIW+W-n^dMTV (22) 

where all constant factors are absorbed in c. It is positive as long as <fi < 2/( (from equa- 
tion 21 and as required for equation 19). Expression 22 must be integrated over N in 
the distribution 16 for all possible peak luminosities L to obtain the differential emission 
measure distribution (DEM), i.e., 

Q(T) = = f S J™ll™ dLo 



\<MT\ J Ll \d\nT\ dL 
ck'T 2 ^ L («/0/(6-«-(«-W(i-/J) dLo 

J Li 



C 



(exponent) 



L 2 
Li 

> T2K L (&-2/C)/(&-^)-(q-2/?)/(i-/?)+i 2 ^ 



where d = ck' is again a numerical constant that is > if (3 < 1 (which is fulfilled; also note 
that (3 7^ 1 is required; see equation 16). The denominator "(exponent)" is identical to the 
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exponent of L which we assume is negative so that large but rare flares do not dominate 
the average emission measure distribution. This holds for reasonable choices of b, (3, and 
4> for a given a but needs to be checked in individual cases. Then, the upper integration 
limit can be set to L 2 = oo as the corresponding term converges to zero, and the complete 
expression 23 remains positive. For the lower limit L\, there are two cases, i) For a given 
T, there is a smallest flare that reaches this temperature at its peak, i.e., T = To, and this 
flare then has a peak luminosity given by equation 13, i.e., L = afT^~^. ii) Since a > 2, a 
lower energy cut-off of the flare energy distribution is required. Flares with energies below 
this limit -E^m are assumed to be unimportant for the heating. The minimum radiated 
flare energy corresponds to a minimum peak luminosity L min = E min /T and was reported in 
Table 4. We therefore require 

Lx(T) = max(a/T fe -*,L mm ). (24) 

Consequently, there are two temperature regimes, depending on whether the first or second 
expression applies. From equation 23, we obtain 

Q oc T 2 ^ ,afT b ^<L min (25) 

K r -(6-*)(a-2/9)/(l-/J)+26-* , a/T ^ > ^ (26) 

The EM distribution thus rises with a slope of 2/£ from low temperatures toward a turnover, 
and then drops again with a power-law exponent that depends on the flare energy distribution 
index a. The turnover occurs where afT b ~* = L min . Note that for C = 0.3-0.7 (Reale et al. 
1997), Q oc T 3 ~ 7 which is much steeper than the prediction from a quasi-static loop model 
but agrees well with measurements of the DEM of the active K star e Eri (Laming et al. 1996) 
and of the late G star £ Boo A (Drake & Kashyap 2001). The DEM thus provides another 
important diagnostic for flare-heated coronae: i) The high-T slope of the distribution can 
be used to determine a; ii) the low-T slope gives information on the heating time-scale of 
flares during their decay; and iii) the peak indicates the turnover energy from a steep to a 
shallow flare energy power-law. - We briefly discuss two cases, for which we use a = 2.2. 

i) For r independent of flare energy, j3 — 0. Our simulations required lowest energies 
corresponding to peak luminosities as low as a few times 10 25 erg s _1 ; we assume L m i n = 
10 25 — 10 26 erg s _1 . The cooling losses per unit emission measure are only weakly dependent on 
temperature in the X-ray range (0.1-10 keV), namely A « 2x 10~ 23 erg cm 3 s" 1 (determined 
from a MEKAL model in XSPEC using our best-fit abundance for Fe). The peak EM required 
for such a flare is thus 5 x 10 47 — 5 x 10 48 cm -3 , which, from the Feldman et al. relation, 
corresponds to a temperature of about 9—16 MK, with considerable scatter. This is indeed 
the temperature regime of peak EM in the emission measure distribution. It is, however, 
unlikely that the flare energy distribution is limited by an abrupt cut-off at low energies. 
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Rather, the power-law index is likely to become smaller at lower energies so that smaller 
flares become less relevant for heating. We therefore set the turnover temperature arbitrarily 
at ~7 MK, in agreement with the observational DEM. The high-T part of the DEM is thus 
determined mainly by plasma between 10—40 MK (see Figure 10). The Feldman et al. 
relation is best fitted by b = 6 ± 1 in that regime, and the cooling function (as calculated in 
XSPEC for an Fe abundance of 0.54 times solar photospheric) has an average slope of = 0.3 
in that range. We propose to set ( = 0.5 — 1, in the light of the above discussion of published 
values. We thus find Q oc T 3±1 for the low-T part of the emission measure distribution, 
and Q oc j^ 2 -")-^ 1 -") = T -o.26+o.36 = T -o.84±o.2o for the high _ T part This model DEM is 

shown in Fig. 10 (dashed) together with a quiescent DEM (dotted) and a DEM from all data 
(solid), both derived from LECS and MECS (using the polynomial and the regularization 
method as described in Kaastra et al. 1996a). Clearly, including the visible flares fills in more 
emission measure at temperatures between 20-50 MK, as expected. This explicitly confirms 
that the high-T DEM is largely related to flaring emission. The fact that the quiescent DEM 
in Figure 10 still contains plasma at ~40 MK is likely to be due to the inclusion of many 
small flare peaks not cut out for this fit. The reconstructed DEM has limited quality given 
the restricted resolution of the BeppoSAX detectors, but the comparison with the theoretical 
DEMs is suggestive except for the presence of a cool peak at 0.15 keV. 

ii) We previously discussed the possibility that (3 = 0.25. Again, we set = 0.3 and 
( = 0.5 — 1. Then, the low-T slope remains the same, 3 ± 1, and the high-T DEM varies like 
Q oc T-°- 2676+0 - 38 = T- L22±0 - 27 , i.e., quite similar to case (i). 

7. Conclusions 

We have investigated the role of statistical flares in coronal heating of magnetically ac- 
tive stars. Long observations of AD Leo were obtained in order to maximize flare statistics. 
Flares have been suspected to play an important role in coronal energy release and subsequent 
impulsive heating of chromospheric material to high temperatures. Chromospheric evapo- 
ration induced by chromospheric overpressure lifts the hot plasma into the corona where it 
fills closed magnetic loops. Since (solar) flares are always related not only to an increase 
in emission measure but to a significant increase in the average plasma temperature, they 
are natural candidates to heat perhaps all of the detected coronal plasma. Recent progress 
in solar physics (Krucker & Benz 1998; Aschwanden et al. 2000; Parnell & Jupp 2000) has 
added new momentum to this hypothesis. 

Active (but quiescent) stellar coronae exhibit a number of features unknown to the non- 
flaring Sun but suspiciously reminiscent of solar (or stellar) flares: i) Very high temperatures 
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Fig. 10. — Emission measure distribution of the quiescent emission of AD Leo, derived from 
the combined BeppoSAX LECS and MECS data. The bins are equidistant in dlogT. The 
emission measures are averages derived from a polynomial fit and a regularization method 
(for details, see Kaastra et al. 1996b). The dashed lines illustrate the approximate slopes of 
the DEM expected from a superposition of flares (arbitrarily adjusted in EM peak ). 
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up to 2—3 keV, similar to temperatures of large solar flares; ii) accompanying, strong non- 
thermal gyrosynchrotron radio emission attributed to relativistic electrons accelerated in the 
initial phase of the flare energy release (Giidel 1994 and references therein) ; iii) high densities 
(> 10 10 cm~ 3 ) reminiscent of (solar) flare densities (Giidel et al. 2001a,b); iv) "anomalous" 
elemental abundances tentatively ascribed to the action of flares, perhaps analogs to solar 
Ne- and S-rich flares (Brinkman et al. 2001; Drake et al. 2001); v) and finally, the presence of 
a large number of strong flares, where the rate of detected flares correlates with the quiescent 
emission level (Audard et al. 2000). 

We have studied the distribution of EUV and X-ray flares in energy, seeking power laws 
of the form dN/dE = kE~ a where k is the normalization of the distribution and a determines 
the steepness of the distribution. We have applied two methods, one based on the count 
rate distribution of binned data, and the second related to the distribution of arrival-time 
differences of the original photon lists (only for EWE). Despite the fundamentally different 
approaches, the results of both methods are in excellent agreement for the EUVE data and 
indicate a = (2.1 —2.3) ±0.1. Simultaneous X-ray observations obtained with the BeppoSAX 
LECS were treated with the first method. The results again overlap with the EUVE results, 
namely a = 2.4 ± 0.2. Only the MECS data show somewhat shallower distributions, with 
a = (2.0 — 2.2) ± 0.2 which can be explained by the harder sensitivity range of the MECS 
detector, and detection bias in terms of flare temperatures. Our results are compatible with 
the findings of Audard et al. (2000) who applied a flare identification algorithm to explicitly 
record flares and to measure their energies. They are further supported by the findings of 
Kashyap et al. (2002) who study further active stars with the EUVE DS. At first sight, 
our a values support a model in which the complete coronae are heated by a statistical 
distribution of flares, involving flares with energies down to a few times 10 29 erg (of radiated 
energy). Also, a model EM distribution based on the superposition of flares of different peak 
temperatures is compatible with the observed EM distribution. However, before we can 
conclude that flares play an important role in the coronal heating process, we should keep 
in mind the following caveats that stellar observations of the present quality are invariably 
subject to: 

(A) All EUV and X-ray observations refer to the radiated energy from the hot plasma. 
There is considerable ignorance of other energy partitions that also contribute to the flare 
energy budget (Wu et al. 1986): Kinetic energy of the upstreaming plasma; potential energy 
of the lifted plasma; energy in waves; energy in accelerated particles; and energy released 
at longer wavelengths. To interpret our results physically, we adopt the following working 
hypothesis: i) The energy initially released in energetic particles is largely thermalized in 
the chromospheric evaporation process, ii) The remaining energy (kinetic and potential) 
is eventually thermalized (e.g., when material drops back to the chromosphere), iii) All 
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thermal energy is eventually radiated away during the cooling processes. We emphasize that 
this applies also to all energy that is conducted from the coronal loops downwards. This 
energy is radiated by the chromosphere, iv) The fraction of the radiative energy released in 
the X-ray range is similar for all flares. 

While points i)-iii) are supported by observations (Dennis 1985) and by numeric simu- 
lations (e.g., Nagai & Emslie 1984; Antonucci et al. 1987), point iv) is relatively difficult to 
assess. Clearly, a considerable part of the energy is lost at UV wavelengths not accessible to 
our observations. The tendency of the flare temperature to increase with overall flare energy 
(Feldman et al. 1995; Aschwanden 1999) would suggest that smaller flares lose a larger frac- 
tion of their energy outside the EUV/X-ray regime. This is, however, of little relevance for 
us: All flares considered here are quite large, with probable peak temperatures (according 
to the Feldman et al. relation) exceeding 20 MK; even the quiescent emission shows its 
peak EM between 5-10 MK (Table 3). Hudson (1991) reports that approximately 2/3 of the 
total radiant energy of a solar flare are emitted in soft X-rays. The total, long-term average 
energy loss in optical U band flares is linearly correlated with the average X-ray losses in 
active stars (Doyle & Butler 1985). If we missed a population of very small flares with very 
low temperatures (e.g., comparable to microflares in the Sun), then they would simply add 
to the flare distribution on the low-energy side, i.e., our distributions would become steeper 
still, making small flares even more crucial for the total energy release. Further, the mea- 
sured average thermal energy input ( "heating rate" ) during a solar flare scales linearly with 
the radiative loss rate in X-rays (Aschwanden et al. 2000). 

Finally, from the phenomenological point of view adopted in the present study, the 
relation between losses in the X-rays/EUV and those at longer wavelengths can be ignored 
altogether if we keep with our goal of modeling the observed coronal emission. The latter is 
clearly dominated by X-ray/EUV radiation. If we successfully explain the total X-ray/EUV 
emission by the radiation of a statistical ensemble of flares, then this simply implies that 
there is no significant additional coronal component. Although there may be additional 
energy release at lower temperatures (i.e., at chromospheric levels), this becomes irrelevant 
for the question of coronal heating. 

(B) The power-law distribution of flares may change spatially on the star. Stellar 
observations unavoidably treat the corona as an average structure. The recent solar results 
with a > 2 were obtained in regions of the quiet Sun while most larger flares occur in active 
regions. It may, however, be interesting to mention that an observation of a single solar 
X-ray bright point resulted in a power-law index similar to those obtained from the whole 
Sun (Shimojo & Shibata 1999). 

(C) The power-law distributions may also vary in time. Bai (1993) and Bromund, 
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McTiernan, & Kane (1995) find a 154 d periodicity in which a changes by ~ 0.2 — 0.4 in 
solar data. There may also be a dependence on the overall magnetic activity level that varies 
with the (cyclic or irregular) "magnetic activity cycle". This latter conjecture is, however, 
not supported by the solar studies of Feldman et al. (1997) and Lu & Hamilton (1991). We 
note that we found different values for a depending on whether or not the early part of the 
DS observation (containing a large flare) was included. 

(D) Although we have used a rather long observing time series (27.3 days of coverage 
with EUVE, referring to segment I-IV), some chance coincidence, like the very large flare 
at the beginning of the observation, may introduce considerable systematic bias. We have 
investigated the role of this large flare on the result for a and found indeed that its selective 
inclusion/exclusion can shift the optimum value by Aa ss 0.1. 

(E) There may be high-energy cut-offs ("roll-overs"; Kucera et al. 1997) related to 
the maximum energy that can be liberated in stellar active regions. In a limited set of 
observations with a limited dynamic range (ratio of strongest to weakest detected flares, also 
depending on the noise level), the deficit of flares close to the high-energy cut-off (because of 
their small occurrence rate) can induce a steepening of a power law. A consequent shallower 
continuation of the distribution toward flares below our detection limit would contribute less 
energy than estimated with our single power-law approach. The present data do not allow 
us to judge on the presence or absence of high-energy cut-offs. The good representation 
by single-power-law flare distributions does presently not argue for their presence. The 
continuation of the power law from detected flares to energies below the detection threshold 
has been explicitly assumed in our energy estimates, and this is no different from any previous 
(solar or stellar) study. For sufficiently small energies, the large number of small flares 
involved begin to overlap in time (the "confusion limit", already evident in our light curves). 
They can no longer be measured individually unless spatial resolution is available. 

(F) Appreciable non-flare contributions to the EUV/X-ray variability are possible (e.g., 
evolution of non-flaring active regions, newly emerged magnetic regions, rotational modula- 
tion of active regions). They would normally add to the low-level variability and may thus 
tend to steepen the count rate distributions. 

Despite these caveats, some of which will be difficult or impossible to avoid in future 
observations, we presently see no compelling argument against our basic finding, namely, that 
flares statistically contribute an important part to the overall coronal radiative losses, and 
that they are therefore good candidates for the coronal heating process per se in magnetically 
active stars. 

Our values for a are very similar to those measured for microflares in the Sun (Krucker & 
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Benz 1998; Parnell & Jupp 2000) despite the 6 orders of magnitude larger energies involved. 
This factor in energy may partly reflect the level of magnetic activity. If so, then the role 
played by microflares in the Sun is played by the much larger flares relevant here in active 
stellar coronae. 

We find independent support for flare heating of active stellar coronae in their coro- 
nal emission measure distribution. By statistically co-adding flaring emission measures by 
weighting them with the dwell time at a given temperature, we derived an analytical ex- 
pression for the differential emission measure distribution. The DEM is characterized by 
a steeply rising low-temperature part and a falling high-temperature part. The slopes and 
the turnover temperature are in principle determined by the flare energy power-law index 
a, the low-energy break in the distribution, and the flare heating parameter ( during the 
flare decay. Previously published DEMs of active stars (e.g., Laming et al. 1996; Drake & 
Kashyap 2001) and the DEM derived here show characteristic shapes that are compatible 
with our expression but are not supported by quasi-static loop models. We suggest that 
the coronal DEMs directly reflect the operation of heating and cooling mechanisms during 
stochastic flares (Gudel 1997; Giidel et al. 1997). 

We conclude this presentation by emphasizing two observational circumstances: i) It 
may be pivotal in which energy range relevant for coronal losses the observations are made. 
Observations that exclusively record the harder part of soft X-rays selectively favor detec- 
tions of large flares and suppress the relevance of low-energy flares (due to the Feldman 
et al. relation; see also discussion in Porter et al. 1995). As is to be expected from the 
flare-heating hypothesis, the quiescent emission is comparatively soft and is therefore also 
underrepresented in hard observations. We have marginally found this effect in our MECS 
observations. One may wonder whether a similar effect exists for non-thermal hard X-rays 
often used for solar flare energy statistics. If they are generated overproportionally in larger 
flares (as suggested by the "Big Flare Syndrome", Kahler 1982, but also by recent obser- 
vations finding that microflares are radio-poor, i.e., relatively weak in the production of 
accelerated particles, Krucker & Benz 2000) then the statistical distributions may be biased 
toward too low a. We have selected the energy range in which the dominant fraction of 
the coronal flare energy is radiated. Also, the efficiency (ratio between observed count rate 
and incident flux) of the DS and the LECS detectors shows only a weak temperature de- 
pendence, i.e., the observations are equally sensitive to plasma over a wide range of relevant 
temperatures, ii) The power-law distribution may depend on the flare energy range con- 
sidered. There are indications in solar observations to this effect, and we can safely state 
that the power laws found here cannot be extrapolated to arbitrary energies: There must 
be a low-energy break (possibly changing to a shallower distribution) in order to confine the 
total radiated power, and there must be a high-energy limit, corresponding to the largest 
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physically possible flares. 
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